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Abstract 

We construct a spinfoam model for Yang-Mills theory coupled to quantum gravity in three dimensional 
riemannian spacetime. We define the partition function of the coupled system as a power series in go^G 
that can be evaluated order by order using grasping rules and the recoupling theory. With respect to 
previous attempts in the literature, this model assigns the dynamical variables of gravity and Yang-Mills 
theory to the same simplices of the spinfoam, and it thus provides transition amplitudes for the spin 
network states of the canonical theory. For SU(2) Yang-Mills theory we show explicitly that the partition 
function has a semiclassical limit given by the Regge discretization of the classical Yang-Mills action. 
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1 Introduction 

The large spin limit has proven to be a useful tool to study the semiclassical physics of the spinfoam 
formalism for loop quantum gravity [T] . Recent developments include the right spacetime dependence of the 
free graviton propagator [5] . On a fixed spinfoam, this limit is related to the emergence of Regge calculus [3] , a 
discrete approximation to General Relativity (GR). This result is very desirable, as all physically interesting 
quantum field theories have useful discrete approximations [J. In particular, this limit is transparent in 
the Ponzano- Regge (PR) model for three dimensional riemannian quantum gravity [5]: the weights of the 
spinfoam partition function are, in the large spin limit, exponentials of the Regge action. Similar results 
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are expected in four dimensfons [Sl[7]. If this is the correct way of studying scmiclassical physics, then it 
should be appHed not only to the partition function, but also to expectation values and to the coupling with 
matter. In [S] it was indeed shown that expectation values of geometric quantities in the PR model tend 
exactly to their classical values in this limit. In this paper we continue along this line of investigation, and 
we construct a coupled spinfoam model of GR and Yang-Mills (YM) theory in 3d riemannian spacetime that 
we show has the right semiclassical limit. 

The relevance of the result concerns also the consistence of the proposed coupling of matter fields to 
spinfoam quantum gravity. In fact, the literature on the argument is somewhat limited!^ Models have been 
constructed specifically for YM fields [13], but they suffer from the following problem: the discretization 
procedure at the root of the quantization treats the gravitational and YM connections in a different way. 
The advantage of doing so is that the geometric information of the quantum gravity sector becomes directly 
available to define the metric dependence of the YM action to be quantized. The drawback is that it leads 
to a spinfoam model where the dynamical variables for GR and YM are attached inhomogeneously to the 
triangulation discretizing the spacetime manifold, thus the model does not provide transition amplitudes for 
the canonical theory: the data carried by the spin network states of the canonical theory would not match 
the data carried by a 2d slice of the spinfoam. 

In this paper we show how one can consistently couple YM theory discretizing both the gravitational 
and YM holonomies in the same way. The key to do so is the use of the generating functional techniques 
introduced in [13]. The same procedure has been recently applied to couple fermionic fields to PR [T5]. As 
a consequence of this procedure, the metric dependence of the YM sector is given by a quantum operator 
whose action can be evaluated using the recoupling theory for SU(2). We study in detail this operator and 
show that the model has the right semiclassical limit but also non trivial quantum corrections. The resulting 
model describes spinfoams with labels attached homogeneously, thus it provides transition amplitudes for 
the canonical theory. 

2 The classical theory 
2.1 Action for GR 

We represent spacetime by a 3d differentiable manifold M, and we use the triad field as fundamental variable 
for GR, as this is more suitable to the spinfoam formalism. The triad is an SU(2)-valued 1-form e^(a;) related 
to the metric by 

gf,u{x) = el^{x)el{x)6ij. (1) 

Geometrically, the triad is a homomorphism between the tangent bundle T{M) and a principal bundle 
P(M, SU(2)) whose fibres take values in the SU(2) groupH Over the bundle P(M, SU(2)) we define a 
connection 1-form ojj/{x), with curvature F{u!). 
We use the following first order action for GR, 

Scnieiix), iol-^x)] = ^ J^^ Tr e A F{lo). (2) 

The trace Tr is over the algebra indices. The variation with respect to the connection gives the torsion-free 
condition, d^je = 0, as an equation of motion. Therefore to is the spin connection. The variation with 
respect to the triad gives F{uj) ~ 0. The theory is invariant under diffeomorphisms acting on the (greek 
letters) spacetime indices, and SU(2) transformations acting on the (latin letters) algebra indices. Under 
these transformations, all solutions of the equations of motion are equivalent, thus the theory has no local 
degrees of freedom [TTj. Indeed 3d GR is the topological BF theory with SU(2) as structure group [T5] . 

We work with dimensionless coordinates, so that lo and F are dimensionless and the triad field has 
dimensions of a length. We use units c = 1, but we keep both h and G. The physical dimensions of the 

^On the other hand, there are very promising results on the coupling of particles [9l 1101 . See 1111 for the coupling of matter 
in the canonical formalism, and 1121 in the group field theory approach. 

^In principle, the triad could take values in SO(3), which is the proper structure group of pure 3d riemannian GR. The group 
is extended to SU(2) ~ Z2 X SO(3), in order to allow the coupling with fermions. Note also that while the bundle P{M, SU{2)) 
is necessarily trivial, there are two isomorphism classes of the bundle P{M, SO{3)), characterized by their Stiefel- Whitney 
number, only one of which is isomorphic to the tangent bundle of M. For a discussion of a possible role of the topology of the 
bundle, as well as the constants h and G, see 1161 . 
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gravitational constant G depend on the spacetime dimension n. In SI units, we have [G] ~ kg^^m"^'^. In 
3d [G] — kg~^ and thus ^ has the correct physical dimensions of an action, namely kg-m, in units c = 1. 
We define the 3d Planck length as £p = SnhG. 

2.2 Actions for YM 

Consider YM theory with gauge group G = SU(A^). We call A° the YM connection. The index a is in the 
algebra g of Q. The YM curvature is F^^ — d^^A^ — duA''^ + f°''"^A^^A'^, where /"''^ are the gauge group 
structure constants. Using the metric as the variable for the gravitational field, the Yang-Mills action on 
n-dimensional curved space is 

Sym h.. ,Al] = ^ J d-x^g^P .9^^ F;, F^, . (3) 

Here g = detg^j^. The YM connection is dimensionless, thus the coupling constant gQ has dimensions 
of an inverse action times {n — 4)/2 length dimensions (in order for the action ^ to have the correct 
physical dimensions for every spacetime dimension n). In the 3d case, we have [h~^go~'^] = [hG]. Conse- 
quentely we introduce a constant with dimensions of a length, called "YM length", ^ym = ft^^ffo"^- The 
dimensionality of the YM coupling constant in 3d is well known, and it is for instance an indication of 
the super-renormalizability of this theory. As it will become clear below, the only relevant constant of the 
coupled system GR plus YM is the adimensional ratio ip/lyu- 

To describe the coupled system with a unique set of variables, we recast this action in terms of the triad 
field (HI). This can be done using twice the relation 

e^^^detel{x) = ^_suK e'^^" e^ix) e^x) (x) (4) 

to write ([3]) as 

SyM[ei{x),A';,{x)] = -^ [ d^xe-'e'^'^'e'^^eleiF^^F:,. (5) 

Notice that no inverse triad appears in the action ([5]), at the price of having to deal with the inverse 
determinant e^^. 

The definition = 6Sym[s, gives the equations of motion, 

F^^(c.) - hlFic) = 16^GT^ d^ei - d^e^ + lj' j,ei - uj' j,e'^ = 0, D^{A) (eF^^A)) = 0. 

As it is well known, the presence of matter does not add degrees of freedom to GR [T7]. Namely the 
equations of motion can be solved without local degrees of freedom in the gravitational sector. Local degrees 
of freedom only come from the Yang-Mills sector of the theory. Geometrically, we are working with a bundle 
P(M, SU(2), Q), where the fibres are independent. This follows from the fact that the YM field A° does not 
couple to the connection ujj/ . 

The coupled system of gravity and YM theory can be described using the action S[e,uj,A] given by ([2]) 
plus where ([2|) describes gravity alone, while ([5]) describes both the YM field and the interaction with 
gravity. To achieve a more symmetric description, we consider an alternative formulation of YM theory, 
which makes use of a first order action: YM theory can be written as a deformation of BF theory [19j 

i J d^x e^'^^BlF^^iA) -^Jd'x ^ggP'^B^B:. (6) 

The auxiliary B field has dimensions of an action. The classical equivalence with ^ can be easily checked 
by looking at the following equation of motion, 

0=^ = ^"^-"^°'^^- 

and using ef^P" e"^^ g^^ = ggf^^g"^^. 

At the quantum level, the formal equivalence can be straightforwardly verified performing the gaussian 
integration over the B field in the partition function for On the other hand, the equivalence of the 
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perturbative quantization of © and ([5]) is rather subtle: particular care has to be taken in the gauge fixing 
procedure |201 121j . These are relevant issues that will however not concern us in the following, where we 
focus only on defining the full spinfoam partition function and studying its semiclassical limit. 

Introducing the Hodge star defined by *i? — ■^^e^upB^dx'^dx'' the interacting terms can be written as 
B A *B, and the coupled action 

S[el{x),Luj/ix),B-^ix),Al{x)] = -L. f Tr e A i^H + / Tr S A - ^ / TvBA^B. (7) 

With abuse of notation, Tr denotes the trace over both algebraic indices / and a. As for ([5]), the metric- 
dependent interacting term can be written in terms of triads only using ([¥]) twice, 

BA*B = ^g^'^ B; Bl = i e'^ e^^^ e^^^ e{ ei B^ B^. (8) 

This is the action which we use in this paper to study quantum gravity coupled to YM theory: it has two 
free pieces, which correspond to topological theories, and an interacting piece with local degrees of freedom. 



3 Discretisation 

In the following, we are going to quantise the theory using functional integrals regularized by means of a 
discretisation of the spacetime manifold^ To do so, we introduce an oriented triangulation A. This is an 
abstract cellular complex made out of points p, segments s, triangles t and tetrahedra r, together with an 
operation that allows the identification of the n — 1 dimensional boundary of an n dimensional object in A. 
To this abstract complex is associated its dual A*, with a one-to-one correspondence between n-dimensional 
objects of A and 3 — n dimensional objects of A*. The 2-skeleton J of A* is a 2-complex made of called 
vertices w, edges e and faces /. 

The coupled action ([7]) is invariant under diffeomorphisms of M . Introducing a fixed triangulation breaks 
this symmetry, thus a procedure to restore it will be necessary in the final theory. A possible way of doing so 
is to use the group field theory formalism |23| . This formalism restores the original diffeomorphism invariance 
by producing in general a sum over triangulations weighted by the corresponding spinfoam amplitude for 
each fixed triangulation. We postpone the constuction of such a model for further work, and we focus here 
on the preliminary step of defining the model on a fixed triangulation. 

On the abstract cellular complex A. we introduce the algebra variables 

Xi e sw(2), b: e fl, 

associated to the segments, and the group variables 

gt e SU(2), Ut e 
associated to the triangles. Furthermore, we define the following quantities: 

e'^^Vllgu e'^^^rllU,. (9) 

t3s t3s 

Here V means that the product over all the triangles which share the same segment s is oriented. The 
orientation is induced from the orientation of s. 

Consider the classical (discrete) theory defined by the action 

S[Xl gt,B^, Ut] - ft^ Tr [X, Z,] + Tr [B, W,] - hXSBB[Xi, B^]. (10) 

s s 

The interaction term is given by 

hXSBB[Xi,B:] = i.n-V^^^lM, (11) 

^Notice that a lattice discretization of the action Q on fiat spacetime has been considered in |22 |. 
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where we introduced the shorthand notations 

(here g^'i^sss jg ^]^g completely antisymmetric tensor for the three segments in a given p, with the convention 
gi23 _ 2 for a right-handed triple) and 

V.[X] ^UKXixiX^, (13) 

where 123 is a right-handed triple of segments for each point p. 

We kept the constants out of the definition of 5bb for a better clarity in later computations. The coupling 
constant between gravity and YM is 

The action PT]]) defines a theory on A, with variables X^ , gt, ,Ut- In the rest of this section, we will 
show that it is an approximation to the continuum theory with action ([7]). Recall that the latter is defined on 
M, with variables e^(a;), ujj/ {x), B^{x), A'^{x). Consider an embedding l : A ^ AI, which allows us to think 
of A as a cellular decomposition of M. Using the embedding, we have ~ da;''. For each point p in each 
tetrahedron r, we call £'^{p) the three vectors tangent to the three segments s belonging to p{t). We can 
think of i'^ip) as a 3 by 3 matrix, and we choose coordinates such that detf^(p) = ^e''^''^''^ef,^pi'^J''^JP^ = 1. 
The embedding i pushes forward to an embedding for the dual triangulation A* , and wc analogously have 

^ dx^ for each edge in A*. 

We make the following four identifications. 

1. The triad field: 

(x) ^ Xl ^ ^el{x)£^, ~ ^ ^ e^(a;)dx^ (15) 

For each segment s in A, we clearly have ll = = £p^XfXf . Therefore the variables Xs 

represent the segment lengths. We will see below that in the quantum theory they are obscrvables 
with a discrete spectrum. This is the sense in which spinfoams describe a quantum geometry, with 
discrete properties for the geometrical observables. 

2. The connection field: 

uj'/{x)^gt = e<''^ ^e^^^'\ (16) 

On the continuum, the connection is geometrically interpreted as an infinitesimal parallel transport. 
On a discrete setting, this should be properly taken over by a discrete minimal parallel transport. The 
group element defined above realizes this property. 

3. The field B: 

B^ix) ^ B1 ^ ^B;{x)i^ - ^-jBl{x)dx^^. (17) 
The geometric interpretation of this auxiliary field is to provide lengths in the fibre space. 

4. The YM connection: 

^;:(x)^C/tEEe^2^= ^e^=^°. (18) 

As for the gravitational connection, the YM connection is discrctised as a group element realizing a 
finite, though minimal, parallel transport along the edges of the dual triangulation. 

We claim that for each configuration of the continuous fields we can find an embedding l such that the 
difference between ([7]) and (fTO]) is arbitrarily small, if the variables in pO)) are interpreted as (fT5|) - p8)) . 

To see it, we start applying the Stokes' theorem to the connection around a closed face / (for non-abelian 
groups, the theorem holds up to corrections in the area of the face): 



E 



oj^t>: ^ I oj^ F{lo). 

df Jf 
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As a consequence of the equality above, the definitions and ^TE\i endow the group variables defined in 
([5]) with the interpretation of holononiies, 



These group elements appear in the first two pieces of pI7|) . When the embedding is sufficently refined, and 
the coordinate areas consequently small, we can expand the group elements around the algebra, so that 
Zg ~ Jj: F{u!) and the first two pieces of ([TU)) reduce to the two BF terms of ([7]). 

Consider now the interacting term Let us choose a point p in a tetrahedron, with 123 a right-handed 
triple of segments. Using (fT5)) and det£'^{p) = 1, we immediately have 



(19) 



where in the last step wc used the fact that e is the volume of the parallelepiped, and there are 3! tetrahedra 
in a parallelepiped. Symmetrizing the expression above over the four points of t we get (jl3p . Analogously, 
using Uni), ini) and det^^(p) = 1 we have 



EpSiS2S3 Ft"- ^ pPf"^ pi pi 



siep(T) 



(20) 



Squaring this and then symmetrizing over the four points, we obtain (fl3 

We summarize the discretisation of the dynamical variables in the following box. 



GR 


GR 




YM 


YM 


smooth variable 


discrete variable 




smooth variable 


discrete variable 


triad 


e su{2) 




auxiliary field 5° 




connection u;^'^ ^ 


ge e SU(2) 




connection A", 


UeeG 



Let us add two remarks before proceeding with the quantization. 

• We have chosen to discrctise both the gravitational and the YM connection on the dual triangulations. 
An alternative procedure, already considered in the literature [13j . is to discretise the YM connection 
on the segments of the original triangulation A, namely as as opposed to A"^l'^ used in (fT8|) . Since 
the gravity sector provides the lengths Is of A, this procedure has the advantage that the geometric 
information for the YM sector is directly described. The resulting spinfoam model presents however 
an awkward feature: the degrees of freedom of YM and GR are attached to different simplices. This 
obstructs a clear connection with the canonical formalism. The discretization proposed here, on the 
other hand, leads to a homogeneous description, as it will become clear below. 

• The YM holonomy Ut enters only the topological part of YM: the interaction term does not contain 
the group elements representing the gravitational and YM holonomies, but only the algebra elements 
Xl and B". This is a consequence of using the action ([7]) to describe the coupled system that will be 
useful to define the spinfoam partition function. 



4 Spinfoams and quantisation of BF theory 

The quantum theory of ([5|) can be constructed from the spinfoam partition function 



t JSU(2) \ J5u(2) 



(21) 



This quantity can be evaluated using the harmonic analysis of SU(2) (for details, see for instance [Tl[9]), and 
one obtains 

=End™j.n{6.7}, (22) 
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where the sum is over all possible assignments of half-integers j to the segments of A. The half-integers, or 
spins, label the irreducible representations of SU(2). The quantity dimj = 2j + 1 is the dimension of the 
representation j. Finally, a symbol is associated with each tetrahedron t of A. The {6j} symbol is the 
key object of the recoupling theory of SU(2), and it depends only on the six js attached to the segments 
of the tetrahedron. It is defined in terms of the Wigner 3m symbols (namely - up to normalization - 
Clebsch-Gordan coefficients) as 





j2 


.73 


H 


I 


3b 


J6 , 





Jl J2 -^S W Jl J5 -^6 W 3i J5 -^3 ] [ Ji n J6 

mi m2 J \ nii mg J \ 7714 TO5 J \ 7714 TO2 mg 



(23) 



where the sums over the repeated m^'s are understood. (j22p defines the Ponzano-Regge model for riemannian 
3d quantum gravity in the absence of matter. 

To quantize the coupled system, notice that in ([7]) there is also a second BF action, corresponding to 
the topological term of the YM action ([6]). For SU(2) YM theory, the quantization of this action is the 
same as in the PR model, only the interpretation of the variables differ. For Q = SU(A^) with N > 2 the 
quantization goes along the same lines described above, but the different algebraic properties of the group 
influence the final expression. In particular in (|23p . one uses a peculiarity of SU(2). namely the uniqueness 
of the Clebsch-Gordan decomposition. The tensor product of two SU(2) irreps is reducible into a direct 
sum of irreps, by means of the formula ji (g) j2 = bi ~ ^2! ffi ■ ■ • ffi ji + J2- The peculiarity here is that each 
irrep in the RHS appears only onceQ This uniqueness is lost for SU(-/V) with N > 2. Consider for instance 
SU(3), and let us use the Cartan highest weight notation to label the irreps; the tensor product of the two 
fundamental representations (1, 0) ® (0, 1), contains twice the adjoint irrep (1, 1). As a consequence, the six 
labels on the segments are not enough to completely characterize a gauge invariant state. The additional 
numbers required are called intertwiners, and they are carried by the 3m symbols. In the SU(3) example 
above, they distinguish between the two adjoints (1,1), which have the same irrep labels. Recalling that 
combinatorially the 3m symbols are associated to triangles, it should be clear that this additional data can 
be used to label the triangles of A. Let us indicate by {it} a collection of intertwiners associated with the 
triangles A. Call p the label for the irreps of SU(A^); the partition function generalizes for Q = SU(-/V) 
to 

Zbf[G]^ lldimp,'[[Arips,it)- (24) 

{Ps,it} T 

Only the labels for the segments and triangles belonging to r enter the tetrahedron amplitude Ar{ps,it)- 
When N = 2, the tetrahedron amplitude reduces to the {6j} symbol described above. In the abelian 
case U(l), the tetrahedron amplitude is trivial; calling the irrep labels, we have only the combinatorial 
condition '^(X^set'^s) ^'^^ each triangle, which is the lattice equivalent of the Gauss law. 

4.1 Expectation values 

Geometric observables in the PR model are gauge-invariant functions of the variables . To compute their 
expectation values it is convenient to introduce a generating functional, 

Z[J] = n / dg^Hf dX, ,^12.T.x.iz.+j.)_ (25) 

t JSV{2) , J,u{2) 

This can be evaluated as described in [Mj, to giv^ 




m = E ridimj. n n i f • (26) 

{js} s ^ \ser 



^The uniqueness of the SU(2) decomposition only concernes the 3d case: for spacetime dimension n > 3 there are generahzed 
3n — j symbols appearing, which carry additional n — 3 quantum numbers. Therefore, in dimension higher than 3, we have 
intertwiners also for SU{2). 

^ To be more precise, the generating functional II25I I is evaluated using a slightly different discretization procedure, in which 
the Jf- variables are not associated with the segments (which are dual to faces), but with the wedges, 2-surfaces introduced in 
| 24l intersecting dual faces and triangles. For the applications of I I25I I in the rest of the paper, the wedge variables would give 
the same results as the X^, thus for simplicity here we do not use them. 
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The quantities entering this expression are defined as follows [T3|. P{J) is a function that relates the Lebesgue 
measure on and the Haar measure on SU(2). Parametrizing as usual the algebra elements as J = ijju, 
with u a unit vector on the 2-sphcrc. we have | J| = "0 and 



P(J) 

The symbols with source insertions is given by 



2 . |J| 



(27) 




Ji 
nil 



32 
1712 



J3 
TO3 



Jl 

ni 



Je 
me 



mil 



J5 



J3 



D^i) fe'^ML'^-'"^ fe'^^ll?^^'^^ (e-'^)D'-^^'^ fe'^-*)£>(^=^ (t 

^miriiV^ )-^m2n2\^ J-^m^nsy^ J-^m4n4\'^ j-^m^n^V 



Ji 
714 



J2 

(e'^«), 



J6 



(28) 



where the Z?'s are representation matrices. The J's are attached to the segments of A; they are the sources 
of the quantum excitations js- If Js = for all segments, reduces to the expression for the {6j} symbol 

given above and P{0) = 1, thus Z[J] = Zp^. 



Z[J] 



(29) 



Using the generating functional the expectation value of a gauge- invariant observable can be 

written as 

(I') = n / dg, n / dX, mi] ^= = * [-^4: 

To evaluate expressions of this type we need to know the action of the algebra derivatives on the measure 
term (|27|) and on the {6j} with sources ([28|) . To evaluate the first action, notice that odd derivatives of 
P{J) always vanish in J = 0, and 



a(2 



9 J^i . . . dJ'^ 



-P{J) 



{-If 



,gll-l2„ 



(30) 



j=o 22"(2n+l) 

where S^^ -'^" is the fully symmetrized tensor, namely S" for ti = 1, (S'-^S^^ + S^^S'^^ + 5'^5''^)/Z for 
n = 2, and so on. These derivatives only produce constant shifts, independent of the irrep label, and thus 
negligible at leading order in the large spin limit. The second action can be evaluated noticing that acting 
on a group element in the representation j, we have 

SJ^ ^ ' j=o 

where T^'^' is the I-th generator in the representation j . By inspecting (j28p . we see that a derivative acting 
on Js attaches to the segment s an algebra generator in the irrep js labeling the segment. This action is 
called "grasping", and it is common in the spinfoam literature (sec for instance [25j [Mj [8] ) . In particular, 
to study (fTT|) we need the action of the double grasping, which enters twice ([T2l) . and the triple grasping, 
which enters the cubic term (jl3p . 

Let us introduce the following diagrammatic notation. 



(31) 



SJL SJL 



Sl 



S2 



SJ 



1 * S S 5 



SJ( SJ^ SJ^ 



(32) 



Notice that to shorten the notation, in the second diagram wc included the sum, which wc recall is over 
the four points in the tetrahedron, and 123 is a right-handed triple for each point. Concerning the double 
grasping, below we will need the cases when si and S2 coincide or share a point. In these two elementary 
cases the graspings give respectively (see [5] jl 



J 



C\j) 




c''in)-c'in)~c'U2) 




(33) 



® The shift i in the double grasping on the same segment was omitted in |8], where we were mainly interested in studying 
the leading order of the semiclassical limit, to which it does not contribute. However notice that this term authomatically 
leads to the Ponzano-Regge ansatz l'^ = (j + i)^ = j{j -I- 1) + ^. I thank Etera Livine for pointing out to me this additional 
contribution. 
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where C^{j) = j{j + 1) is the SU(2) Casimir. We do not report the rather lengthy resuh for the triple 
grasping, for which we refer to [S] . Notice that the only contribution from ([3D]) to these graspings is the shift 
of 1/4 to the double grasping on the same segment. For the expectation values on a single tetrahedron r, 
we introduce also the following notation, 



Si 



S2 





(34) 



We proceed in the same way for the YM sector, introducing a source 77° for the B'^ field and evaluating the 
generating functional as before, keeping in mind that for Q = SU(A'^) the partition function for BF depends 
also on additional quantum numbers labeling triangles, the intertwiners. Defining the SU(-/V) equivalent of 
([28]) (where the intertwiners label the 3m symbols), we have 



{Ps,'t} 



To distinguish the GR graspings from the YM ones, we picture the latter with curly lines, such as 

Sl 

<5 (5 i 



As above, we have the elementary graspings 



P2 



S2 



p pi' * p''^ pi' 

Notice that now the non diagonal grasping depends also on the intertwiner label i, and / depends on the 
gauge group considered. For Q = SU(2) there is no intertwiner and f{p) ~ ^ \C^{p^) — C^(pi) — C^(p2)] as 
before. 

Taking the product of ((26l) and pS]) we define the generating functional for the two BF theories. 



P2 



» PS 



(35) 



(36) 



(37) 




(38) 



Using this generating functional we can compute the expectation values of functions of and B^. In 
particular of (|12p which gives the interaction between GR and YM. With the double graspings for GR and 
YM defined above the expectation value of (fT2|) can be computed from the action of the following operator. 



Cr 



SJL SJt{ SJi 5Ji Srj-^ 



(39) 



where the diagrammatic notation implicitly includes the sums, so the lines represent all the possible graspings 
within the tetrahedron. However the action of this operator on ([35]) is ambiguous, as it contains the product 
of two GR graspings: for non abelian groups the graspings do not commute, so to properly define this 
operator an ordering prescription is needed. We do so in the next section, motivated by recovering the right 
semiclassical limit. In section [6] below we come back to this generating functional and construct the spinfoam 
model for the coupled system. 



5 The semiclassical limit 

The spinfoam model that we are going to construct is defined on the fixed triangulation of spacetime 
introduced in section [3j Therefore we do not expect its semiclassical limit to reproduce the continuum 
action ([2]), but a suitable discretization of it. For the Ponzano-Regge model ([^^ this is well understood: in 
the limit when all the spins are large the discrete approximation to GR described by Regge calculus emerges. 
In this section we briefly recall the way this happens, and use it to fix the ordering prescription needed to 
deal with products of graspings. 
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5.1 Regge calculus 



It is an old result [3] that on a triangulation A, made of flat tetrahedra, the curvature is distributional and 
concentrated on the segments. All the geometric information is encoded in the segment lengths is- The 
action can be written in terms of the segment lengths and their deficit angles e^, 

S[ls\ = Y,^ses{ls)- (40) 

The sum here is over all segments in A. The dependence of the deficit angles on the segment lengths can be 
found as follows. First, the deficit angles are defined in terms of the dihedral angles as eg = 27r — X^ras ^1 ■ 
Then, the dihedral angles of a tetrahedron can be expressed in terms of the segments, through the well 
known formula 




where Ai and A2 are the areas of the two triangles ti and ±2 sharing the segment s, as shown in the above 
figure. Finally the areas A = V{2) and the volumes V = V(3) can be expressed in terms of the segment 
lengths as determinants of the Cayley matrix, 



(n) 



(-!)"+! 

— — — det C(„\, 
2"(n!)2 



C(„) = 



V 



1 





1 



1 \ 



/2 



(42) 



/ 



A non trivial result of Regge calculus is that the deficit angles are the conjugate variables to the segment 
lengths, 

Einstein's equations are read from dS/Slg = on the bulk. From ((43|) . we see that they imply a flat 
triangulation, with zero deficit angles everywhere in the bulk. 

To see how this formalism for classical GR emerges from the PR model, we consider the Regge action on 
a single tetrahedron, 

SK[is]=Y,'^-^siQ- (44) 

The key fact is that exponentials of this action dominate the (homogcneout0) large spin limit of the {6j} 
symbol entering ([22]) [3 [J 

cos(^R(j,) + f) 



(45) 



\/l27ry(j,) 

where V{js) and S^{js) are respectively the classical volume (|^ and the Regge action of the tetrahedron 
with segment lengths given by ^ = £p (j + i). The factor j does not change the equations of motion, and 
the presence of the cosine as opposed to a single exponential is due to the fact that the PR model sums over 
both orientations of the tetrahedron, and the action ((44|) is odd under change of orientation (see discussions 
in 

If taking this limit is the correct way to study semiclassical physics, then also the geometric quantities 
that one can evaluate in the model should reduce to their classical expressions. Indeed, in it was shown 
that the double and triple graspings reduce to classical scalar products and volumes. For the double grasping 
on two segments i and j sharing a point, we have the asymptotics 




Vl27ry(j, 



''Namely we rescale the half-integers entering the {6j} symbol as js = Nks and take the A'^ 
notation, we do not write explicitly this rescaling but keep using the variables js- 

*This feature will also complicate the analysis of the coupled model as we discuss below. 



(46) 



00 limit. To keep a simple 
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The scalar product here is defined in terms of the segment lengths in the natural way: if denotes the 



third segment in the triangle defined by ig^ and £s2i we have 2£si ■ = ^'t^ ~ ^si ~ ^'to- ^'^^ the triple 
grasping, we have identically zero if the three segments are coplanar, and in any other case we get 




(47) 



Notice that the volume is purely imaginary, and it has a different phase in the argument of the cosine. As 
described in [8] and anticipated in [l^ , this is due to the fact that the PR model sums over both orientations 
of spacetime, and the definition (|13p of volume changes sign under change of orientationlf] Consequently, 
real results with the correct phase can be obtained looking at even powers. 

As the graspings do not commute with each other, ordering ambiguites arise in defining the products 
of graspings. As discussed in [8], we can choose a "temporal ordering" T prescribing the graspings to be 
performed one after the other. Under this prescription we have 



T{ 




is J 



cos(gR(j.) + f) 
n/12^F(Js) 



(48) 



T{ 



} 




V{js 



3S (SrUs) 



if k is odd, 



if k is even. 



(49) 



This is the way Regge calculus emerges from spinfoams in the large spin limit: the partition function 
reduces to exponential of the Regge action, and the expectation values of the geometry to their classical 
counterparts. This means that the theory has a well defined semiclassical limit. Below we are going to 
construct a spinfoam model for the coupled system of GR and YM. If this framework is indeed robust, then 
the coupled system should have a semiclassical limit, given by large irrep labels, in which gravity is described 
by Regge calculus and YM theory by a discretization a la Regge, which we construct below. 



5.2 Regge discretization of YM theory 

We discuss here how ([SJ can be written in the Regge approach introduced above. First of all, we consider 
the topological BF term. To keep things explicit, let us focus on the special case G = SU(2). By analogy 
with (|45|) . we know that the large p limit of p4|) will be dominated by the cosine of the discrete action 
S = Bsil's{Bs) where Bs = fi{ps + ^) represents lengths in the fibre space. The "dihedral angles" 
ips measure the curvature of the fibre space, and the equations of motion impose it fiat. This is just a 
straighforward discrete version of the BF action. More interesting is the Regge description of the term ([5]) 
which makes YM dynamical and carries the interaction between GR and YM. We discretize this by a sum 
over flat tetrahedra as 

^ J Sx^gg^'^BlBl^ig^Y^Vrg^'BlBl, (50) 

where we have replaced the canonical volume form (Px^fg by (six times) the tetrahedron volume Vr- To 
express (|50|) in terms of Regge variables, we need the flat inverse metric as a function of the segment lengths. 
To obtain this expression, we we start from the formula for the volume V of the tetrahedron, 

(3! Vf = .9 = ^ e^" -''" £'^"••'^"5^1., • • • (51) 

Deriving this, we obtain 

£ " i ^""'"^ ■ ■ " ^^^^ 

^This could be avoided introducing a modulus in the definition of the volume, as is done in the canonical approach |25| . 
However, this would prevent a definition of the functional derivatives in terms of graspings. 
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Expressing the derivative in terms of the squared edge lengths = gfj.i^l'^i^, we read 

where the sum is over the six segments in the tetrahedron. 
Using ([55)) and (fT7|) , we can write ([50)1 as 

5Rint[4, B,] = 3! h'go' E E a^T ^s- (54) 

To compute the volume derivatives, it is convenient to introduce a double index notation, where a segment 
is identified by its vertices, s = ij with i,j = 1 . . .4. Let us also introduce the vectors normal to the 
triangle obtained removing the point i from the tetrahedron. We choose their orientation so that they point 
inward. These vectors satisfy = AAf, fii ■ ftj = AAiAj cosffy , where 9ij are the dihedral angles defined 
in (|1T|). Then, deriving ioi n — 3 and using (pij) we obtain 



dV 1 n, 



(55) 



del 72 V ■ 

In conclusion, we can write the interaction term on a Regge triangulation as 

^Rlnt[4, B,] = 1 h'g„' E E ^ Bty (56) 

We have defined the ordering prescription to deal with the products of graspings. and described the 
Regge description of the classical coupled system. Below we will show that this is exactly what emerges in 
the semiclassical limit of the quantum theory. 



6 Spinfoam model of the coupled system 

We now come back to the construction of the spinfoam model for the coupled system. To define the quantum 
theory for the coupled system, we evaluate the partition function 

Firstly, we regularize this formal expression by means of the discretization procedure discussed above, ob- 
taining 



^GRYM = IT / dgtT] f dUtT] f dX,Y\ f dB, 

f'JsU{2) fJQ f'J5u{2) s-^S 



J^^[TrX,Z,+Ti-B,Ws]-i\SBB[X,B] 



(58) 



with A ^BB defined in . We see explicitly from the expression above the advantage of the symmetric 
form of the action obtained using the first order formalism for YM theory. In fact, the interaction term 
'S'bb [X, B] is a gauge invariant function of the algebra variables only, and does not contain gravitational nor 
YM holonomies. We can thus use the generating functional ([38]) and with the substitutions i-^ ~^sT''^ 

-i-^rl^ write 

^GRYM - e-'^^^-[-'^'-'^lz[J, r;] 

J=?)=0 

We can interpret the equation above in the following sense: Z[J,rj\ gives a colouring of the triangulation, 
with labels for SU(2) and Q irreps on the edges, and labels for Q intertwiners on the triangles. This colouring 
provides the kinematical arena for the dynamics of YM and the interaction between YM theory and gravity, 
which is realized by the action of the local operator e*'*''^^^. 
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To define Sbb[ 



the inverse volume (V[X]) 
define this operator as 



], we use the T-ordering introduced above, but we also need a prescription for 
appearing in pip . Motivated by recovering the right semiclassical limit, we 





6 ' 











dt 



OO 

k=0 



i-tl 
kl 



6J 



2k+l 



(59) 



where t is an auxiliary variable with no physical meaning and the powers of the triple grasping are T- 
ordered as in The advantage of this definition is that, using (P^ . we have in the large spin limit the 

desired semiclassical behaviour of the inverse volume, 



SJ 



1 COS (SrUs) + |7r) 



V{js) ^UttVUs) 

With these definitions, we can evaluate the partition function as a power series in A: 



where 



^1 



BB 



^GRYM 



5 J' 6t] 



E 



s 



BB 



.5 .5 

6 J' Sr] 



J=T]=0 



(60) 



(61) 



dt 



OO 

E 

k=Q 



kl 



Vr 



SJ 



2k+l 



Cr 



The expression (|6ip provides us with a formula for the partition function (|57p in the spinfoam formalism. It 
is well defined, order by order in A. The ordering prescription for the operators is as in (j48p and (j49p . Sbb 
is local because the graspings only connect segments belonging to the same tetrahedron. Notice that this 
partition function provides a quantization of YM theory only when the full series is considered. As discussed 
in [21j . the order- by-order equivalence requires additional gauge- fixing terms, which change the topological 
nature of the zeroth order and thus the spinfoam procedure used here. 

The zeroth order of the coupled partition function (|6ip corresponds to the product of two independent 
BF partition functions, 

= Z[J,r,] = Zbf[SU{2)] ZbfIG]. (62) 

J=r)=0 

It is a sum over all the labelings of all the segments of the triangulation by irreps of SU(2) and Q, where 
the SU(2) labels describe the geometry of the triangulation, as in the pure gravity picture, and the Q labels 
describe the degrees of freedom of YM theory. We report in the following table the homogeneity of the 
labeling, which has been achieved discretizing both GR and YM connections on the dual triangulation: 



simplex 


GR labels 


YM labels 


GR amplitude 


YM amplitude 


Interaction 


point 












segment 


3 


P 


dim j 


dim p 




triangle 




i 








tetrahedron 






m 


Ar{p,i) 


SBB[js,Ps,it] 



This structure allows us to identify 2d slices of the spinfoam with a spin network state for the canonical 
quantization of the coupled system. The links of such a spin network would be coloured by the labels for the 
irreps of GR's SU(2) and YM's SU(iV), and on the nodes by the intertwiners of SU(A^): \s) I7, j;, p;, i„). 
In the next section we study the metric dependence of the interaction term S'BB[is7 Ps-, it]- 

^^li we compare it with the Turaev-Viro model 1261 . we see that t behaves like the contribution, from a single tetrahedron, 
to an imaginary cosmological constant, A = it. 
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6.1 Effective action 



The first order term in the expansion (j6ip is Z'^' = —iX{SBB), thus it amounts to computing the effective 
action for the interacting term ([8]). For the SU(2) graspings we use the diagrammatic notation p2l) . and 
we T-order them as in (|48p . (j49p . The double grasping on Q has no ordering ambiguities (it appears only 
linearly at this order, and it commutes with the other graspings). Introducing the diagrammatic notation 
as in (|32p . we can write 




The key diagram to evaluate is 



{Cr) = T{ 



}: 




(64) 



This diagram has 4 x 36 = 144 contributions, coming from all the possible choices of graspings in a given 
point, times the four points. Each contribution can be evaluated using grasping rules and recoupling theory 
as in [8]. Because there are only double graspings entering this expression, the evaluation is rather simple 
and we do not report the details here, but only the asymptotics. 

Let us distinguish two types of terms, when the YM grasping is diagonal, namely S3 = ^3, and when is 
not diagonal, namely S3 7^ ^3. Consider first the diagonal case. For fixed S3 — ij, there are 4 contributions 
from p = i, and four from p = j. To fix ideas, let us choose S3 = 12. The YM grasping immediately gives 
C'^(pi2) from (|37p . As for the GR grasping, we can implement the antisymmctrization from the e tensors of 
([M)) in the algebraic indeces, and using the large spin limit of the relevant graspings from p = \ gives 



' yJ T'^l 



COS 



44 - (^13 • ^14)' 



(65) 



The term in square brackets can be immediately recognized as 4^2^. Analogously, the graspings from p = 2 
give 4:Ai^. There are in total 48 contributions of this type. 

The second case is when the YM grasping is non diagonal. Let us fix S3 = 12, ^3 = 13. Choosing 
for simphcity SU(2), the YM grasping gives ^[C^(p23) — C^{pi2) — ^^(^13)]. The GR graspings have four 
contributions all coming from p = 1, times 2 from the symmetric choice S3 = 13, ts = 12, giving in total the 
eight contributions 



4 ( -^^13 X(l 



cos [S^jjs) + f) 



?12-^14)(43-^14)-(^12-^13)4 



(66) 



If we recall the definition of the dihedral angles 9s in terms of angles 4>ij between the segment vectors (defined 
by (.ij ■ £ik = iijiik cos^jfc), 



sm sm (pjk cos ffik = cos cpij cos (pjk - cos (pik , 



(67) 



the term in square brackets in ([66]) reads n2 ■ n^. Notice that the same contribution comes also from the 
choice S3 = 24, ^3 = 34, for which the YM grasping gives ^[C^(p23) — C'^{p24) — C^(p34)]. There are in total 
96 contributions of this type. 

Adding up all the contributions we get 



1 



' 4 n2 • ^3 



2 ^C2(P23) - C\p,2) - C2(pi3)) + 2 (C2(P23) - C2(p24) - ^^(^34) 



(68) 
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Collecting all the Casimirs and using the fact that X]i=i ' 
Therefore has the following large spin behaviour, 



adds up to -8 J^i-, ■ C^(py ). 



cos {Snifs 



(69) 



Next, we look at the other diagrams entering ([55)) . Thanks to the T-ordering, the powers of the triple 
grasping factorize in the large spin limit and one obtains simply (|60p . Putting everything together and 
approximating the sums over the irrep labels with integrals, we then have 



-iA(S'BB) 



.90 • rij ^2 



cos (SrUs) 



Itt) cos(S'r(ps 



(70) 



In the limit in which C'^{pij) ^ BfjEf^, the term in square brackets is exactly ■^S'Rint given in 156]) . This is 
our key result: the large spin behaviour of the spinfoam model (j6ip is dominated by the classical dynamics 
of GR coupled to YM described a la Regge. 

6.2 Tetrahedral asymptotics of the coupled partition function 

Let us now go back to the full partition function (|6ip . and consider a single tetrahedron. In this simple case, 
we can proceed as above to show the general result for powers of the operator ^bb , 



{Sr{Ps 



i cos (S'R(js) + Itt) if n is odd. 



cos (S'R,(js) + f ) if n is even. 



(71) 



Using (|7ip we can sum up the series (j6ip and evaluate the coupled partition function in the large spin limit 

cos(5R(p,0 + f) 



J Yl '^^s '^Ps n ^"^y s dimps 



^l2TiV{]s)^l2nV{ps) 



E 



(*5RI„t)^" + ^ 



^ (2n 



1)! 



E 

n=0 



RInt ) 

(2^ 



\2n 



COS \S^{js) + - 



W dj, dp. n dim^ dimp, eos(gR(,.)-5RR,(j.,P.) + f) cos (5r(p..) + f ) 

We see from this result that the asymptotics of the tetrahedral partition function are dominated by linear 
combinations of exponentials of the classical discrete actions described in section O Notice that in principle 
one would expect a single exponential of the Regge version of ([7]) from a correct semiclassical limit, namely 



Yet the presence of the somewhat akward linear combination in (j72p can be completely understood if we 
recall that the spinfoam quantization includes a sum over both orientation of the (triangulated) spacetime 
manifold. In fact, let us look at (jT]): under change of orientation (e^ i— *■ — e^) both the first and third terms 
change sign, but not the second. This accounts for the first cosine term in ([72]) . 



The second cosine comes naturally when the same argument is applied also to the fibre manifold: under 
change of orientation of the fibres {B^ '-^ —B^) only the second term in ([7]) changes sign, whereas the first 
and third do not. 
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So the tetrahedral partition function (|72p is the correct semiclassical Umit of a coupled quantum theory 
with the feature of summing over both orientations of the spacetime and fibre manifold. However it is 
exactly the presence of this sum over orientations that makes problematic the interpretation of the partition 
function on a generic triangulation, whose semiclassical limit we do not report here. In the PR model this is 
not a problem: the model is topological invariant, thus the partition function on a generic triangulation can 
be reduced to the one on a single tetrahedron. The model presented here on the other hand is manifestly 
non topological invariant, thus this argument does not apply. We leave this question open, stressing that 
it is an issue of the PR model in its own, and not of the coupling with YM performed here. Furthermore, 
notice that the problem can be circumvented by introducing oriented boundary states in the computation 
of physical correlations, as it is done in the graviton calculations of p]. 

7 Conclusions 

In this paper we defined a spinfoam model for riemannian 3d quantum gravity coupled to Yang-Mills theory. 
The partition function of this model is given in ((6T|) . The construction makes use of the generating functional 
technique, and of grasping rules and SU(2) recoupling thery. With respect to previous attempts in the 
literature, the model has the advantage of discretizing the gravitational and YM connections in the same 
way. This leads to a homogeneous labeling of the triangulation by the GR and YM variables, as reported in 
the table in section [51 and thus this model provides transition amplitudes for spin networks of the canonical 
theory. Furthermore, for the particular choice of SU(2) YM theory, we were able to show explicitly that the 
model has the correct semiclassical limit, given by the discretization a la Regge of the classical YM action. 
This result supports both the consistency of coupling matter fields to gravity in the way proposed here, and 
the utility of the large spin limit to study semiclassical physics. 

Let us add a few comments and perspectives on further developments. One thing to stress is that as 
usual, the quantization procedure presents many ambiguities. In particular, here we chose a definite way 
of discretizing the dynamical variables and the actions, and we gave a particular ordering prescription for 
the products of non commuting operators. Other choices are possible. This one was preferred because it 
easily reproduced the correct semiclassical limit. Once the correct limit is established, the key question 
becomes computing the quantum gravity corrections to it. Already in the simple model presented here we 
see that this is a formidable task due to the large number of contributions that the next to leading order 
of the asymptotics receives. Yet this crucial question certainly deserves further work. Among the intriguing 
effects to look for there are modified dispersion relations: even though the discrete structure of loop quantum 
gravity does not break Lorentz invariance per se [27], it is commonly recognized that possible deformations 
of this symmetry is a promising sector to look for observable quantum gravity effects [28| . Remarkably in 
[TU] it was shown that the coupling of point particles to the Ponzano-Regge model gives a well defined non 
commutative effective action for the particles, thus resulting in a precise Lorentz deformation with modified 
dispersion relations. It would be extremely interesting if analogue effects arise in the model here presented 
A potential difference lies in the fact that here the matter field is added to the gravitational action, whereas 
in the approach of |TU] the particles are described using the gauge degrees of freedom of the gravitational 
sector. However the situation could be similar to what happens in quantum field theory in the temporal 
gauge, for both YM theory [29] and linearized gravity [30], where static matter sources (whose coupling is 
lost by the choice of gauge) turn out to be described by gauge degrees of freedom. 

At the end of the previous section we mentioned how an oriented boundary state could fix the orientation 
of the coupled partition function. Indeed, finding a suitable boundary state would have far more important 
applications: it would allow to follow the proposal developed in [5] to compute scattering amplitudes for the 
coupled system described here. This would be a key development to understand low energy physics in the 
context of spinfoam gravity. From this point of view, the emergence of Regge calculus in the large spin limit 
is very promising, as this property is at the basis of the results in [2] , and can be used to construct a suitable 
boundary state, as discussed in [31j . 

Another important question for future work is how to restore the full diffeomorphism invariance of the 
continuum theory, here broken by the choice of a fixed triangulation. Being the model non topological, the 
symmetry should be restored by including a sum over the triangulations, for instance along the lines of the 
group field theory approach |23| . Notice that restoring the right degrees of freedom can only be expected 

Notice that the quantum particle of 3d YM theory is massless and spinless, thus the results of 1101 might not apply 
straighforwardly. 
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for the fully rcsummcd partition function. This is because as pointed out in section 6, only the full partition 
function can be expected to be equivalent to YM theory. This important aspect will be studied elsewhere. 

Finally, it is sometimes argued that the spinfoam quantization of SU(2) BF theory simply amounts to 
mapping the B field to its eigenvalue j. This map has been considered in the literature, especially in 4d, to 
define the spinfoam quantization of additional terms to the BF action, may they be constraints reducing BF 
to GR, or matter fields. We would like to stress that particular care should be taken in using this map: the 
more consistent procedure used in this paper, where the additional terms act as genuine quantum operators, 
shows that such a map would fail to capture the full quantum theory. In particular, one would not see 
the ordering ambiguities and would fail to reproduce the full richness of the quantum corrections to the 
semiclassical leading order. 
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